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The theory of the thermal hydrodynamic fluctuations in bounded systems has been 
developed only for non-superfiuid liquids on the base of the linearized Navier-Stokes equa- 
tions [1]. However, the study of effect of confined geometry on correlation functions of 
superfiuid liquid is important for explanation of experiments in which capillaries, pores, 
dense suspensions and other objects are probed by radiation [2]. 

In this paper we present the correlation theory of the equilibrium thermal hydrody- 
namic fluctuations of compressible superfiuid helium confined by the spherical cavity of 
radius R. The hydrodynamic fluctuations in the superfiuid helium are described on the 
base of linearized two-liquids Khalatnikov's equations [3] with spontaneous heat fluxes 
and strains. For mass flow 6 J and relative velocity 6w = Vn — Vg of normal-fluid and 
superfiuid components these equations have form 

^+V-5J = 0, p^+p.sV.5w-^A5T = -iv-g, 

^+V5p-^A5J-i(ir^+C2)V(V-5J)-^A5^+^(pCi-Jr^-C2)V(V-M = -/, 

ot p p 6 p p 6 

-^-Pn^+'^Sp-psV6T-CiV{V-6J) + ^i-pC,+p%)V{V-6w) = -pVH, (1) 
ot ot p 

where 
6p = 

Here s is the entropy per unit mass, p is the pressure; k, t], and Ci)C2,C3 are the heat 
conductivity, the "first" viscosity, and the "second" viscosity coefficients, respectively; 
Ps,Pn, and p are the superfiuid, normal-fluid, and total densities. The system (1) is 
written with account of the Onsager symmetry principle for the kinetic coefficients. All 
kinetic coefficients and thermodynamic derivatives in (1) assumed to be independent of 
coordinates. By (1) the fluctuation fields of mass flow, mass density and other variables 
are linear functionals of random forces in the right parts of these equations. 





Assuming that the rigid-cavity surface is motionless we apply zero boundary conditions 

— * 

for the perpendicular component of mass flow 5 J and temperature fluctuations 5T in the 
liquid helium, and the stick boundary conditions for the normal component of velocity [4] 

5Jr{f,t) = 0, Vn{f,t) = 0, ST{f,t) = (r = R). (2) 

We present the mass flow, relative velocity, mass density and the temperature fluctu- 
ations in the form 

SJ{r, ^) = E {j'it)Lx{r^ + Jxit)M,{r-) + {t)Nx{f)) , (3) 

A 

5wif, = E {w^{t)L,ir) + W^f (t)M,(r) + {t)N,{r)) , (4) 

A 

5p{r,t) ^Y.Mt)jn{kr)Y^n{0,^), (5) 

A 

ST{f,t) = Y.Tx{t)jn{kr)Ymn{e,ip), (6) 

A 

where A = {n, Z, m} and 

Lmn{r) = T-^{Ymn{0,v)3n{kr)), 

Mmn{r) = rot{fY^n{0, (fi)jn{kr)), Nmn{r) = ^rotMmnir)- (7) 

ft^A 

Here Yynn{d, ^) are the spherical functions and jn{kxr) are the spherical Bessel functions 
[5]. Note that L^„(r), M^„(r), iV^„(r) form the system of pairwise orthogonal functions. 
Prom the boundary conditions (2) we find out eigenvalues k\: 

jn{Pnl) = 0, kx = Pnl/R, (8) 

Inij'nilni) + ajn(7nO =0, -n,n + l, kx^ Inl/R- (9) 

— * 

We present also the longitudinal and transverse parts of the random viscous force / 
and heat fluxes q in the form of expansion (3) with coefficients and Q\, respectively, 
a — M, N, L, and the random potential if - as in eq.(5) with coefficients Hx. 

It is convenient to use the Fourier representation for random SJ, Sw, Sp, ST 

oo 

x{t) = J x^exp[—iut]duj. (10) 

— oo 

Since the equations for the Fourier- amplitudes of the mass flow J^^ and relative velocity 
Wx^ fluctuations {b — M, N) do not contain components of other fields we immediately 
obtain that 

rM,N 

J.. - PnWx^ - _ . (11) 



Ignoring the products of dissipation coefficients, fluctuation forces and keeping only the 
terms linear in the dissipation coefficients we get 
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+iuj I klpsS 
where 
Dx = -pnP 
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gfj kxUj + iajxUj){LJ + kxUj + iajxUj) 



(12) 



(13) 



and Ui, U2 are the first- and second-sound velocities for the infinite space [3]; aix, Oi2x are 
the attenuation coefficients of the first- and second-sound, respectively. Ignoring terms 
proportional to ps and — l) ^ for the first- and second-sound attenuation 

coefficients we obtain 
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where Cp, c„ are the specific heats at constant pressure and volume, respectively. The 
expansion coefficients for the excess mass density can be expressed in terms of J^. 

From the expression for the rate of energy dissipation in the liquid Hell confined in 
the sphere volume V 



(15) 



Q{t) 
we find that 
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dV ( -Vnf + Ps6w WH - -5T V • q 



Q.^E (--{PsWl + JD^lfl + hxPsW^A^^Hx^ , a^L,M,N, (16) 

a,X \ P / 



where A" = / df\a{r)\'^, a — M,N,L are the normalization factors. 



V 



Using the Langevin equations (11), (12) we construct as in [6] the matrix of the gener- 
alized susceptibilities. And then by FDT find the spectral densities of equilibrium thermal 
fiuctuations for the expansion amplitudes of the hydrodynamic fields. This resulted in 
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(17) 



where G(a;,T) is the expression [6] for the average energy of the quantum oscillator. 

The result differs from the spectral densities for an infinite superfluid by only the 
discrete behaviour of the wave numbers kx- Factorisation of the coefficient Dx in (17) 
allows us to decompose the expressions for the pair correlation function of into the 
simple fractions. And the spectral densities of the Jx as in the case of an infinite medium 
[7] , consist of two contributions 

I .L|2 ^ 2pQ{iJ,T)uj'^{klul - u^) j {aixui - a2AM2)(^^ - kluj) - aixUikl{uj - uj) 
' ' ' T^^{ki{ul - ulY \ [{u - kxu,f + alxul] [{uj + kxu.f + 

{a2\U2 - aixui){uj'^ - kluf) + a2\U2kl{ul - ul) \ ^ 
[{u: - kxU2Y + alxul][{uj + kxU2Y + alxul] /' " 
Now it is easy to construct the spectral densities of thermal fluctuations of the mass 
flow and the mass density hydrodynamic flelds 

< 5J{f,t)5J{f',t') < \-Jx\' M^Alif'), A^M,L,N, (19) 

X,A 

< 5p{f, t)Sp{f', t') < l^^r >u> jn{kxr)jn{kxr')Y^n{d, f)Ymn{d', f'), (20) 

X 

where 

< \Rx\' 4 < ■ (21) 
We have separated the flrst- and second-sound contributions to the mass density - 

mass density correlation function. The correlation functions are presented as the sum of 
the series in different spherical harmonics and roots of the transcendental equations (8), 
(9). As mentioned above, the difference between this case and that of inflnite space is 
only in the discrete nature of the wave numbers kx- Each term in this series generates 
two Mandelstam-Brillouin doublets, which are caused by quasiphonons of two kinds, one 
having the first- and the other one the second-sound speed. Each components of the 
doublets are of Lorentz type, their widths being given by aixUi and a2xU2, and the spectral 
densities of the thermal fiuctuations of mass fiow and mass density can be presented as 
combinations of infinite numbers of such Lorentzians. 
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